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Abstract 

In the barrier resilience problem (introduced Kumar et al, Wireless Networks 2007), we are 
given a collection of regions of the plane, acting as obstacles, and we would like to remove the 
C*~) minimum number of regions so that two fixed points can be connected without crossing any 

region. In this paper, we show that the problem is NP-hard when the regions are fat (even 
when they are axis- aligned rectangles of aspect ratio 1 : (1 + e)). We also show that the problem 
is fixed-parameter tractable (FPT) for such regions. Using our FPT algorithm, we show that if 
the regions are /3-fat and their arrangement has bounded ply A, there is a (1 + e)-approximation 
that runs in 0(2^ A ^^n 7 ) time, where / € 0(^f- log^A/e)). 

1 Introduction 

u The barrier resilience problem asks for the minimum number of spatial regions from a collection T> 

u that need to be removed, such that two given points p and q are in the same connected component 

& of the complement of the union of the remaining regions. This problem was posed originally in 

i— —i 2005 by Kumar et al. |XQ^ ITT] , motivated from sensor networks. In their formulation, the regions 

^1 are unit disks (sensors) in some rectangular strip B C M 2 , where each sensor is able to detect 

movement inside its disk. The question is then how many sensors need to fail before an entity can 
move undetected from one side of the strip to the other (that is, how resilient to failure the sensor 
system is). Kumar et al. present a polynomial time algorithm to compute the resilience in this 
case. They also consider the case where the regions are disks in an annulus, but their approach 
cannot be used in that setting. 

o 

2 1.1 Related Work 

Despite the seemingly small change from a rectangular strip to an annulus, the second problem 
still remains open, even for the case in which regions are unit disks in M 2 . There has been partial 
progress towards settling the question: Bereg and Kirkpatrick [2 J present a factor 5/3 approxima- 
tion algorithm for the unit disk case. Afterwards, Alt et al. [T] and Tseng and Kirkpatrick |14| 
independently show that if the regions are line segments in M 2 , the problem is NP-hard. 

The problem of covering barriers with sensors is very current and has received a lot of attention 
in the sensor network community (e.g. 0111 [8]). ^ n the algorithms community, problems involving 
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Figure 1: (a) A path of minimum resilience (equal to 1). (b) Three paths with minimum thickness (2). 

region intersection graphs are also popular. Gibson et al. [7] study the opposite problem of ours: 
compute the maximum number of disks one can remove such that p and q are still separated. 

1.2 Results 

The NP-hardness constructions in [1] and [2] both rely on the fact that line segments are long 
and skinny. In this paper, we show that it is not the skinnyness of line segments what makes 
the problem hard. We show that when the regions are fat regions [5] in M 2 , the problem is still 
NP-hard, even for the case in which the regions are axis-aligned rectangles of aspect ratio 1 : 1 + e 
and 1 + £ : 1. Our construction suggests that hardness relies on the fact that the regions can cross 
each other, and it does not easily extend to the case where the regions are disks. In Section [3] we 
discuss this in more depth. 

To complement our hardness result, we show in Section [4] that when the regions are fat the 
problem is fixed-parameter tractable on the resilience. In Section [5] we show that, if the system of 
regions has bounded ply (13}, we can use the FPT result to obtain an approximation scheme. In 
particular, the constructive results apply to the original unit disk coverage setting (since unit disks 
are fat regions). Our results rely on tools and techniques from both computational geometry and 
graph theory. 

2 Preliminaries 

We denote with p and q the points that need to be connected, and with T> the set of regions that 
represent the sensors. We say that any collection of objects in the plane are pseudodisks if the 
boundaries of any two of them intersect transversely at most twice. 

We formally define the concepts of resilience and thickness introduced in [2]. The resilience of 
a path it between two points p and q, denoted r(w), is the number of regions of T> intersected by 
7r. Given two points p and q, the resilience of p and q, denoted r(p,q), is the minimum resilience 
over all paths connecting p and q. In other words, the resilience between p and q is the minimum 
number of regions of T> that need to be removed to have a path between p and q that does not 
intersect any region. From now on, we assume that neither p nor q are contained in any region of 
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D. Notice that such a region must always be counted in the minimum resilience paths, hence we 
can ignore them (and update the resilience we obtain accordingly). 

Often it will be useful to refer to the arrangement induced by the regions of T>, which we denote 
by A(T>). Based on this arrangement we define a weighted graph G^m\ as follows: each vertex of 
Gmx>\ corresponds to a cell of A(V). We also add directed edges between any pair of neighboring 
cells. An edge has cost 1 if following the edge we enter an element of T>, or cost if we are leaving 
it instead. 

The thickness of a path tt between p and q, denoted t(ir), equals the number of sensor region 
intersections of tt (possibly counting the same region multiple times). Given two points p and q, 
the thickness of p and q, denoted t(p,q), is the value |sp G (p, q)\ + A(j>), where sp G (p, q) is a 
shortest path in G^x>) from the cell of p to the cell of q, and A(p) equals the number of regions that 
contain p. Note that the thickness between two points can be efficiently computed in polynomial 
time. Also notice that the resilience (or thickness) between only two points only depends on the 
cells to which the points belong to. Hence, we can naturally extend the definitions of thickness to 
encompass two cells of A(V), or a cell of A{T>) and a point p G M 2 . 

Note that thickness and resilience can be different (since entering the same region several times 
has no impact in the resilience, but is counted every time for the thickness), see Fig. [TJ However, 
as we will see later, the thickness (and the associated shortest path) will help us find a path of low 
resilience. 

The following useful lemma follows directly from Lemma 1 in [2], and applies to the case where 
the regions in T> are unit disks. In the following, "well-separated" is in the sense used in [2] (i.e., 
the distance between p and q is at least 2\/3)q 

Lemma 1. Let T> be a set of unit disks, and let S C T> be an optimal solution. If p,q are well- 
separated, then no disk of S appears more than twice. 

Corollary 1. When the regions ofD are unit disks, the thickness between two well- separated points 
is at most twice their resilience. 



3 NP-hardness 

In this section we will show that computing the resilience of certain types of fat regions is NP-hard. 
First we show NP-hardness for general connected regions, and later we extend it to axis- aligned 
rectangles of aspect ratio 1 : 1 + e and 1 + e : 1. We start the section establishing some useful 
graph-theoretical results. 

3.1 Maximum independent set and Euler paths in oddly embeddable graphs 

Let G be a graph, and let p be a point in the plane. Let V be an embedding of G into the plane, 
which behaves properly (vertices go to distinct points, edges are curves that do not meet vertices 
other than their endpoints and do not triple cross), and such that p is not on a vertex or edge of 
the embedding. We say T is an odd embedding around p if it has the following property: every 
cycle of G has odd length if and only if the winding number of the corresponding closed curve in 
the plane in T around p is odd. Fig. [2] shows some examples. We say a graph G is oddly embeddable 
if there exists an odd embedding T for it. 

^Note that the well-separatedness of p and q is only needed to prove a factor 2 rather than 3. Everything still 
works for ill-separated points, at a slight increase of the constants. Our most general theorem statements, for /3-fat 
regions, do not require the points to be well-separated. 
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(a) (b) (c) 

Figure 2: (a) A planar odd embedded graph, (b) A non-planar one. (c) A tripartite graph, oddly embedded 
around p. 

Lemma 2. Every tripartite graph is oddly embeddable. 

Proof. The vertices of a tripartite graph G can be divided into three groups V\ , V2, V3 such that 
there are no internal edges in any of these groups. Now, consider a triangle A around p. We create 
a drawing V where all vertices in V\ are close to one corner of A, the vertices in V2 are close to a 
second corner, and the vertices in V3 are close to the remaining corner. All edges are straight line 



segments. See Fig. 2(c) 



Consider the graph H obtained from G by contracting all vertices in Vi to a single vertex Vi\ H 
is a triangle (or a subgraph of a triangle). Now consider any cycle in G, and project it to H. Since 
there were no edges in G connecting vertices within a group Vi, this does not change the length of 
the cycle, nor does it change the winding number around p. Any two consecutive edges from V{ to 
Vj, and back from Vj to vi, do not influence the parity of the length of the cycle, nor the winding 
number around p, so we can remove them from the cycle. We are left with a cycle of length 3w and 
winding number w or —w, for some integer w. Clearly, 3w is odd if and only if w is odd. Therefore, 
r is an odd embedding of G, as required. □ 

The maximum independent set problem in a graph asks for the largest set of vertices in the 
graph such that no two vertices in the set are connected by an edge. This problem is well-known 
to be NP-hard on general graphs. In fact, it remains NP-hard for tripartite graphs. A simple proof 
is included for completeness, and because we need the argument later. 

Observation 1. Let G = (V,E) be a graph. Let G' be obtained from G by subdividing every edge 
e £ E into an odd number of pieces, by adding an even number m e of new vertices. Let m = ^2 e m e 
be the total number of vertices added. Then G has a maximum independent set I C V if and only 
if G' has a maximum independent set V with \L'\ = \L\ + m/2. 

Proof. For every independent set / C V in G, there is a corresponding independent set I' in G' 
with I L'\ = \L\ + m/2: for every pair of extra vertices on an edge, we can always add one of the 
two to an independent set. Conversely, for every independent set V in G' , there is a corresponding 
independent set I in G with |/| = |/'| — m/2: I cannot use both extra vertices on an edge, 
so if we simply remove all extra vertices we remove at most \E\ elements from /' (clearly, if we 
remove less than m/2 vertices from /' this way, we can remove more vertices until / has the desired 
cardinality) . □ 

Corollary 2. Maximum independent set is also NP-hard on tripartite graphs. 
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(a) (b) (c) 

Figure 3: (a) An oddly embedded graph with four crossings, (b) The crossings are flattened according to 
the parity of their vertices, (c) An odd Euler path from p to the outer face. 



Corollary 3. Maximum independent set is NP-hard on oddly embeddable graphs. 

Given an embedded graph T, we say that a curve in the plane is an odd Euler path if it does 
not go through any vertex of T and it crosses every edge of T an odd number of times. 

Lemma 3. Let p be a point in the plane, and T an oddly embedded graph around p. Then there 
exists an odd Euler path for T that starts at p and ends in the outer face. Moreover, such path can 
be computed in polynomial time. 

Note that in particular, if T is a plane graph, then the lemma follows directly from the conven- 
tional notion of Euler paths in the dual (multi-)graph. 

Proof. First, we insert an even number of extra vertices on every edge of T such that in the 
resulting embedded graph T', every edge crosses at most one other edge. Now we construct an 
Euler path that crosses every edge of V exactly once; note that this path will therefore cross every 
edge of T an odd number of times. Consider a pair of crossing edges and the four vertices concerned. 
For each pair of consecutive vertices (vertices that are not endpoints of the same edge) , find a path 
in the graph that does not go around p (when seen as a cycle, after adding the crossing). 

The parity of the length of this path does not depend on which path we take: if there would 
be an even-length path and an odd-length path between the same two vertices, both of which do 
not go around p, then there would be an odd cycle that does not contain p, which contradicts the 
oddly embeddedness of V . Now, if the path has even length, we identify these two vertices. Note 
that of the four pairs of vertices involved in a crossing (i.e. ignoring the two pairs forming edges in 
r), exactly two pairs will have odd length connecting paths, so effectively we "flatten" the crossing. 
We do this for all crossings, and call the resulting multigraph V" . See Fig. |3(b) (If the two crossing 



edges belong to different connected components of T, there are no paths connecting their vertices; 
in this case we make an arbitrary choice of which vertices to identify.) 

Now r" is planar. Furthermore, by construction, all faces of T" have even length, except the 
one containing p and the outer face. Therefore, the dual multigraph of T" has only two vertices 
of odd degree, and hence has an Euler path between these vertices. Furthermore, this Euler path 
crosses every edge of V exactly once, and therefore every edge of T an odd number of times. Notice 
that the proof is constructive. Moreover, both the transformations and the Euler path can be done 
in polynomial time, hence such path can also be obtained in polynomial time. □ 
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(1)) 



(c) 



Figure 4: Creating regions to follow T and T. 



3.2 Reduction from resilience to maximum independent set 

Lemma 4. Let p be a given point in the plane, and T an oddly embedded graph (not necessarily 
planar) around p. Furthermore, let T be a curve that forms an odd Euler path from p to the outer 
face. Then we can construct a set V of connected regions such that a minimum set of regions from 
T> to remove corresponds exactly to a minimum vertex cover in T. 



Proof. If T is self-intersecting, we can rearrange the pieces between self-intersections to remove 
all self-intersections. Thus we assume that T is a simple path. 

If T crosses any edge of V more than once, we insert an even number of extra vertices on that 
edge such that afterwards, every edge is crossed exactly once. Since we inserted an even number of 
vertices on every edge, by Observation [T] this does not affect the maximum independent set. Let 
r' be the resulting graph. 

Now, for each vertex v in P, we create one region D v in T>. This region consists of the point 
where v is embedded, and the pieces of the edges adjacent to v up to the point where they cross T. 
Fig. |4(b)| shows an example (the regions have been dilated by a small amount for visibility; if the 
embedding T has enough room this does not interfere with the construction). Note that all regions 
are simply connected. 

Finally, we create one more special region W in D that forms a corridor for T. Then W is 
duplicated at least n times to ensure that crossing this "wall" will always be more expensive than 
any other solution. Fig. 4(c)| shows this. 



Now, in order to escape, anyone starting at p must roughly follow T in order to not cross the 
wall. This means that for every edge of P that T passes, one of the regions blocking the path (one 
of the vertices incident to the edge) must be disabled. The smallest number of regions to disable 
to achieve this corresponds to a minimum vertex cover in P. □ 

Recall that minimum vertex cover (or maximum independent set) is NP-hard on oddly embed- 
ded graphs by Corollary [3| As a consequence, we obtain our first hardness result for the barrier 
resilience problem (but note that this also follows from [1] and |14j . although their constructions 
are quite different from ours). 

Theorem 1. The barrier resilience problem for arbitrary connected regions in the plane is NP-hard. 

Using fat regions. We now adapt the approach to also work for a much more restricted class 
of regions: axis-aligned rectangles of sizes 1 x (1 + e) and (1 + e) x 1 for any e > (as long as 
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e depends polynomially on n). For simplicity, we limit T to have maximum degree 3. Maximum 
independent set is still known to be NP-hard in that case ||6J, and making them tripartite does not 
change the maximum degree. 

The idea of the reduction is the following. We start from a sufficiently spacious (but polynomial) 



embedding of T, as illustrated in Fig. |5(a) On each edge we add a large even number of extra 



vertices. Each vertex will be replaced by a rectangle, so every edge in V will become a chain of 



overlapping rectangles, like the green rectangles in Fig. 5(b) Therefore the first phase consists in 
replacing the embedding of V by an equivalent embedding of rectangles. We call these rectangles 
graph rectangles (green in the figures) . Some care must be taken in the placement of graph rectangles 
around degree-3 vertices and in crossings, so that the rest of the construction can be made to work. 
Next, we place wall rectangles (orange in the figures; these consist of many copies of the same 
rectangle) across each graph rectangle. The gaps between adjacent wall rectangles should cover 
the overlapping part of two adjacent graph rectangles, so that a path can pass through them only 
whenever one of the two graph rectangles is removed. Then, we find a curve T from p that goes 
through every gap exactly once (note that T exists, by Lemma [3]). Fig. 5(c) [ illustrates this phase 



of the construction. Finally, we add more wall rectangles around T, to force any potential escape 
path from p that doesn't go through the wall rectangles to be homotopic to T. Fig. 5(d)| shows 



the final set of rectangles. Now, computing an optimal resilience path among this set of rectangles 
would correspond to a maximal independent set in V. 

For the construction to work, there needs to be enough space to place the wall rectangles. It is 
clear that this is possible far away from the graph rectangles, but close to the graph rectangles we 



proceed as follows: first, Fig. 6(a) shows the placement of rectangles along an edge of T. Fig. 6(b) 
shows how to place the rectangles at degree-3 vertices. Crossings are handled as shown in Fig. 6(c) 
These gadgets force some of the gaps in the chain to join each other. But this is no problem if every 
edge has enough rectangles. Also, note that at the center of the construction in Fig. 6(c)| there are 



two overlapping green rectangles, which belong to the two crossing chains. This is the only place 
where we vitally use the fact that the regions are not pseudodisks. 

Lemma 5. Let p be a given point in the plane, and T an oddly embedded graph with maximum 
vertex degree 3 (not necessarily planar) around p. Furthermore, let T be a curve that forms an odd 
Euler path from p to infinity. Then we can construct a set T> of axis-aligned rectangles of aspect 
ratio 1 : (1 + e) such that a minimum set of regions from T> to remove corresponds exactly to a 
minimum vertex cover in V. 

Proof. We first add groups of extra vertices on every edge of V so that we have room to place 
the rectangles, in an even number per edge. Then replace edges by chains as of rectangles as as in 
Fig. [6j and connect the orange (wall) rectangles to force the only optimal path from p to the outer 
face to be along the Euler path T. The path may have to be rerouted locally close to the crossings, 
but since there is a sufficiently large number of crossings with every edge anyway, this is always 
possible. Orange rectangles have to be duplicated sufficiently many times again, to make sure that 
no optimal path will ever cross them. □ 

Theorem 2. The barrier resilience problem for regions that are axis-aligned rectangles of aspect 
ratio 1 : (1 + e) is NP-hard. 

A similar approach can likely be used to show NP-hardness for other specific shapes of regions. 
However, it seems that a vital property is that they need to be able to completely cross each other: 
that is, the regions in T> should not be pseudodisks. Thus, if one were to prove that vertex cover for 
oddly embeddable graphs of bounded degree is NP-hard would also imply that the barrier resilience 
problem for unit disks is also NP-hard. 
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(c) (d) 

Figure 5: (a) A non-planar oddly embedded cubic graph, embedded on a grid, (b) A set of rectangles, 
containing exactly one rectangle for each input vertex, and an even number of rectangles for each input edge. 
Note that crossings can be embedded if sufficiently far apart, (c) Local walls are added to make "tunnels", 
each tunnel contains the overlapping part of two adjacent yellow rectangles. To go through a tunnel, one of 
the two yellow rectangles has to be removed. Then we choose an Euler path from p to the outside, that goes 
through each tunnel exactly once, (d) The final set of rectangles, designed to force any path from p to the 
outside to be homotopically equal to the one we drew. 
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(a) (b) (c) 



Figure 6: Local details of the construction. Note that we use rectangles with a large aspect ratio for 
visibility, but the same constructions can be made with aspect ratio arbitrarily close to 1. (a) Overlapping 
rectangles to create edges (with an even number of extra vertices), (b) A vertex of degree at most 3, which 
is just a single rectangle, (c) A crossing between two chains. 




Fi gure 7: Any minimum vertex cut disconnecting s and t gives a minimum resilience path between p and 
q (that stays in the rectangular region). 



4 Fixed-parameter tractability 

In this section we introduce a fixed-parameter tractable (FPT) algorithm for the case in which D 
is a collection of /3-fat objects. For clarity we first explain the algorithm for the special case of unit 



disks. In Section 4.2 we show how to adapt the solution to fat regions. The parameter of our FPT 
algorithm is the length of the optimal solution, thus our aim is to obtain an algorithm that, given a 
problem instance whose solution has resilience r, runs in 0{2f^n c ) time, for some constant c and 
some polynomial function /. Note that for treating the case of unit disk regions we assume that p 
and q are well-separated, so we can apply Lemma [l} 

First we give a quick overview of the method of Kumar et al. [10J for open belt regions. Their idea 
consists in considering the intersection graph of T> together with two additional artificial vertices 
s,t with some predefined adjacencies. There is a path from the bottom side to the top side of the 
belt if and only if there is no path between s and t in the graph. Hence, computing the resilience 
of the network is equal to finding a minimum cut between s and t (see Fig. [7]). Our approach is to 
find for a long path that passes through p and q, cut open and transform the problem instance into 
one with something similar to an open belt region. We then follow the approach of Kumar et al. 
taking into account that the right and left boundaries of our region correspond to the same point. 
Hence, instead of using a regular vertex cut, we will need to use a vertex multicut |15j . 

Consider the shortest path r between the cells containing p and q in Ga(t>) , let t be the number 
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(a) (b) 

Figure 8: (a) A set of regions T>, and a shortest (minimum thickness) path between p and q. (b) The 
arrangement shaded by thickness from p. 

of traversed disks (recall that we assumed that p is not contained in any region, hence this number 



is exactly the thickness of p and q, see Fig. 8(a) ). We observe that cells with high thickness can be 



ignored when we look for low resilience paths. 

Lemma 6. The minimum resilience path between p and q cannot traverse cells whose thickness to 
p or q is larger than 1.5t. 

Proof. We argue about thickness to p; the argument with respect to q is analogous. Let p be a 
path of minimum resilience between p and q, and let r be the resilience of p. Recall that p does not 
enter a disk more than twice, hence the thickness of p is at most 2r < 2t. Assume, for the sake of 
contradiction, that the thickness of some cell C traversed by p is greater than 1.5i. Let pc be the 
portion of p from C to q. By the triangle inequality, the thickness of pc is less than 0.5i. However, 
by concatenating r and pc we would obtain a path that connects p with C whose thickness is less 
than 1.5i, giving a contradiction. □ 

Let R be the union of the cells of the arrangement that have thickness from p at most 1.5f; we 



also call R the domain of the problem. Fig. 9(a) shows an example. By the previous lemma, cells 



that do not belong to R can be discarded, since they will never belong to a path of resilience r, see 
Fig. [8(b)} 

Note that the number of cells remaining in R might still be quadratic, hence asymptotically 
speaking the instance size has not decreased (the purpose of this pruning will become clear later). 

We extend the shortest path r from q until a point q' on the boundary of the domain. Let r' 
denote the extended path (Fig 9(b) [ ) . 



Lemma 7. There exists a path r' from p to a point q' on the boundary of R via q, whose thickness 
is at most 1.5t. 

Proof. Consider any shortest path tree from p in the dual graph of the reduced domain R, defined 
as the corresponding subgraph of G^x>) • All leaves of the tree correspond to cells on the boundary 
of R, which by definition have thickness at most 1.5t from p. Therefore all other cells in the tree 
lie on a path from p to a boundary cell that has length exactly 1.5t, including q. □ 
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(a) 



(b) 




(c) (d) 

Figure 9: (a) The reduced domain R. (b) A path r' from p to q' via q. (c) After cutting along r', we get 
the domain R' . We add a set S of extra vertices on the boundary of R' , and we now have two copies of q. 
(d) A crossing pattern, consisting of a topological path p (defined by the sequence of points of S it passes), 
and a binary assignment to the disks of T> intersected by r'. 



11 



We "cut open" through t', removing the cut region from our domain. Note that cells that are 



traversed by t' are split by two copies of the same Jordan curve (Fig 9(c) ). After this cut we have 
two paths from p to q. We arbitrarily call them the left and right paths. 

Consider now a minimum resilience path denoted p; let r = r(p) denote its resilience. This 
path can cross r' several times, and it can even coincide with r' in some parts (shared subpaths). 
Although we do not know how and where these crossings occur, we can guess (i.e., try all possibil- 
ities) the topology of p with respect to t'. For each disk that r' passes through, we either remove 
it (at a cost of 1) or we make it an obstacle. That way we explicitly know which of the regions 
traversed by r' could be traversed by p. Additionally, we guess how many times p and r 1 share 
part of their paths (either for a single crossing in one cell, or for a longer shared subpath). For 
each shared subpath, we guess from which cell p arrives and leaves (and if the entry or exit was 
from the left or right path). We call each such configuration a crossing pattern between r' and p. 
Fig. |9(d)| illustrates a crossing pattern. 

Lemma 8. For any problem instance T>, there are at most 2 4r log r +°( r lo s r ) crossing patterns between 
t 1 and p, where r = r{p). 

Proof. First, for all disks in r', we guess whether or not they are also traversed by p. By Lemma[7j 
t' has thickness at most 1.5t, there are at most such many disks (hence up to 2 choices for which 
disks are traversed by p). 

We now bound the number of (maximal) shared subpaths between p and r': recall that p passes 
through exactly r = r{p) disks, and visits each disk at most twice. Hence, there cannot be more 
than 2r shared subpaths. Observe that r' cannot traverse many cells of A(T>): when moving from 
a cell to an adjacent one, we either enter or leave a disk of T>. Since we cannot leave a disk we have 
not entered and r 1 has thickness at most 1.5t, we conclude that at most 3t cells will be traversed 
by r' (other than the starting and ending cells). 

For each shared subpath we must pick two of the cells traversed in r' (as candidates for first and 
last cell in the subpath). By the previous observation there are at most 3t candidates for first and 
last cell (since that is the number of cells traversed by r'). Additionally, for each shared subpath 
we must determine from which side p entered and left the subpath (four options in total). Since 
these choices are independent, in total we have at most 2r x (3t x 3t x 4) 2r = 2r • 36 2r • t options 
for the number of crossing patterns. Combining both bounds and using the fact that t < 2r, we 
obtain: 

2i-5t x 2 r . 3g 2r . £ 4r < 2 5r x 2r • 36 2r • (2r) 4r = 2 9r+1+logr+2rlog36+4T ' logr = 2 4r log r +°( r log r ) 

□ 

Notice that the bound is very loose, since most of the choices will lead to an invalid crossing 
pattern (for example, patterns that forbid passing through a disk, but force p to pass through one 
of its cells). However, the importance of the previous lemma is in the fact that the total number 
of crossing patterns only depends on r (and not on n). 

Our FPT algorithm consists in considering all possible crossing patterns, finding the optimal 
solution for a fixed crossing pattern, and returning the solution of smallest resilience. From now 
on, we assume that a given pattern has been fixed, and we want to obtain the path of smallest 
resilience that satisfies the given pattern. If no path exists, we simply discard it and associate 
infinite resilience to it. 
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(d) (e) (f) 

Figure 10: (a) We may schematically represent W as a circle, since the geometry no longer plays a role. 
Partial paths are shown dashed, (b) The intersection graph of the regions after adding extra vertices for 
boundary pieces between points of S U {p, q}, shown green, (c) The secondary graph H, representing the 
forbidden pairs, (d) A possible solution (the edges to remove are marked in orange), (e) The multicut 
vertices, translated back to the original regions, (f) The final set of regions and a path of resilience zero 
from p to q. 



4.1 Solving the problem for a fixed crossing pattern 

Recall that the crossing pattern gives us information on how to deal with the disks traversed by 
t'. Thus, we remove all cells of the arrangement that contain one or more disks that are forbidden 
to p. Similarly, we remove from T> the disks that p must cross. After this removal, several cells of 
our domain may be merged. 

Since we do not use the geometry, we may represent our domain by a disk W (possibly with 
holes). After the transformation, each remaining region of T> becomes a pseudodisk, and p becomes 
a collection of disjoint partial paths, each of which has its endpoints on the boundary of W (see 



Fig. 10(a)). To solve the subproblem associated with the crossing pattern we must remove the 
minimum number of disks so that all partial paths are feasible. 

We consider the intersection graph Gj between the remaining regions of T>. That is, each vertex 
represents a region of T>, and two vertices are adjacent if and only if their corresponding regions 
intersect. Similarly to [TO], we must augment the graph with boundary vertices. The partial paths 
split the boundary of R into several components. We add a vertex for each component (these 
vertices are called boundary vertices). We connect each such vertex to vertices corresponding to 
pseudodisks that are adjacent to that piece of boundary (Fig. |10(b)[ ). Let G x = (V X ,E X ) be the 
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resulting graph associated to crossing pattern X. Note that no two boundary vertices are adjacent. 

We now create a secondary graph H as follows: the vertices of H are the boundary vertices of 
Gx- We add an edge between two vertices if there is a partial path that separates the vertices in 
Gx (Fig. 10(c)| ) . Two vertices connected by an edge of H are said to form a forbidden pair (each 



partial path that would create the edge is called a witness partial path). We first give a bound on 
the number of forbidden pairs that H can have. 

Lemma 9. Any crossing pattern has at most 2r 2 + r forbidden pairs. 

Proof. First notice that there are only "a few" boundary vertices: since partial paths cannot cross, 
each such path creates a single cut of the domain. This cut introduces a single additional boundary 
vertex (except the first partial path that introduces two vertices). Recall that we can map the 
partial paths to crossings between paths t' and p and, by Lemma [8| these paths can cross at most 
2r times. In particular, there cannot be more than 2r + 1 boundary vertices. Observe that we only 
add a forbidden pair between boundary vertices of Gx, hence in the worst case we have at most 
( 2r + 1 ) = 2r 2 + r forbidden pairs. □ 

The following lemma shows the relationship between the vertex multicut problem and the 
minimum resilience path for a fixed pattern. 

Lemma 10. There are k vertices of Gx whose removal disconnects all forbidden pairs if and only 
if there are k disks in T> whose removal creates a path between p and q that obeys the crossing 
pattern X . 

Proof. Let A' be the regions of A(T>) inside R that are not covered by any disk (after the k disks 
have been removed) , and let R' be their union. Two points are connected if and only if they belong 
to the same region of A'. Hence, there is a path between p and q with the fixed crossing pattern 
if all partial paths are inside R' . The reasoning is analogous to the one used by Kumar et al. [10J. 
If all partial paths are possible, then no forbidden pair can remain, since — by definition — each 
forbidden pair disconnects at least one partial path (the witness path). On the other hand, as soon 
as one forbidden pair remains connected, there must exist at least one partial path (the witness 
path) that crosses the forbidden pair. Thus if a forbidden path is not disconnected, there can be 



no path connecting p and q. Fig. 10(e) illustrates how the resulting set of vertices corresponds to 



the original question. □ 



That is, thanks to Lemma [10[ we can transform the barrier resilience problem to the following 
one: given two graphs G = (V,E), and H = (V,E') on the same vertex set, find a set D C V 
of minimum size so that no pair (u,v) G E' is connected in G \ D. This problem is known as 
the (vertex) multicut problem |15j . Although the problem is known to be NP-hard if \E'\ > 2 
there exist several FPT algorithms on the size of the cut and on the size of the set E' [12] 115] . 
Among others, we distinguish the method of Xiao ([15j. Theorem 5) that solves the vertex multicut 
problem in roughly 0((2k) k+i ^ 2 n 3 ) time, where k is the number of vertices to delete, I = \E'\, and 
n is the number of vertices of G. 

Theorem 3. Let T> be a collection of unit disks in M?, and let p and q be two well-separated points. 
There exists an algorithm to test whether r(p, q) < r, for any value r, and if so, to compute a path 
with that resilience, in 0(2^ r ^n 3 ) time, where f(r) = r 2 logr + o(r 2 logr). 

Proof. Recall that our algorithm considers all possible crossings between p and r'. For any fixed 
crossing pattern X, our algorithm computes Gx-, and all associated forbidden pairs. We then 



execute Xiao's FPT algorithm |15] for solving the vertex multicut problem. By Lemma 10 the 
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Figure 11: A /3-fat region D is contained in a big disk, but contains a smaller disk; in this example, (3 — 2. 



number of removed vertices (plus the number of disks that were forced to be deleted by X) will 
give the minimum resilience associated with X. Note that the disks to be deleted can be retrieved 
as well. 

We now focus on the running time. Computationally speaking, the most expensive part of the 
algorithm is running an instance of the vertex multicut problem for each possible crossing pattern. 
Observe that the parameters k and I of the vertex multicut problem are bounded by functions of r 
as follows: k < r and i < 2r 2 + r (the first claim is direct from the definition of resilience, and the 
second one follows from Lemma j9j). Hence, a single instance of the vertex multicut problem will 
need 0((2r) r +( 2r2+r )/ 2 n 3 ) = 0{2( T +^ r ^ r2+1 - 5r ^n 3 ) = q^ 2 io g r+o(r 2 iogr) n 3) time> By Lemma|8|the 
number of crossing patterns is bounded by 2 4r log r +°(' r lo § r ) . Thus, by multiplying both expressions 
we obtain the bound on the running time, and the theorem is shown. □ 

We remark that the main focus of this research is the fact that an FPT algorithm exists. Hence, 
although the dependency on r is high, we emphasize that the bounds are rather loose. 



4.2 Extension to Fat Regions 

We now generalize the algorithm to similarly-sized /3-fat regions. A region D is /3-fat if there exist 
two concentric disks C and C whose radii differ by at most a factor (3, such that C C D C C' 
(whenever the constant (3 is not important, the region D is simply called fat). Since we need the 
regions to be of similar size, we assume without loss of generality that the radius of C is 1 and the 



radius of C is (3; in this case we will call D a /3-fat unit region. Fig. 11 shows an example of a 2-fat 
region. We first need to extend Lemma [T] 

Lemma 11. Let V be a set of f3-fat unit regions, and let S C V be an optimal solution. In the 
sequence of disks of S found when going from p to q (in an optimal way), no disk of S appears 
more than (2/3 + l) 2 times. 



Proof. Let D be a region in S, and consider its containing disk C , with center c. Clearly, the 
number of times an optimal path visits D is upper-bounded by the number of times it visits C. 
Now, analogously to the original argument by Bereg and Kirkpatrick |2], every time the optimal 
path visits and leaves C , it must do so in order to avoid some other region. This other region must 
intersect D, and since it is /3-(unit)fat, it must contains a unit disk centered at distance at most /3 
from D. Therefore all regions intersecting D have their unit-disks centered at distance at most 2/3 
from c. Moreover, they are totally contained in a disk of radius 2/3 + 1 centered at c. A simple area 
argument shows that at most (2(3 + l) 2 disjoint unit-disks fit into a disk of radius (2/3 + 1). □ 
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Figure 12: Example showing that the 0(/3 2 ) bound is tight. 



Fig. 12 illustrates that a path may indeed visit the same region as many as £l((3 2 ) times, so the 



bound in Lemma 11 is asymptotically tight. 



Corollary 4. When the regions ofT> are /3-fat unit regions, the thickness between two points is at 
most (2/3 + l) 2 times their resilience. 

This change in the upper bound of the thickness in terms of the resilience implies similar changes 
in Lemmas [6j [7j [8] and |9j The following lemmas summarize these changes; they are proved in the 
same way as their counterparts for disks, thus we only sketch the differences with the original proofs 
(if any). 



Lemma 12. The minimum resilience path between p and q cannot traverse cells whose thickness 

)- 



to p or q is larger than (1 + (2/3 + l) 2 " 1 - 



Proof. (Sketch) We use the same reasoning as in the proof of Lemma [6j On the one hand there 
is the minimum thickness path between p and q, whose thickness is t. On the other hand, we also 
have the minimum resilience path p between the same points, whose thickness is at most (2/3 + l) 2 t. 
Assume now that any cell C traversed by p has thickness k(2/3 + l) 2 i from p, for some < k < 1. 
The alternative path goes from p to C, via q, and has thickness (1 — fc)(2/3 + l) 2 t + t. The bound 
we need is obtained for the value of k that makes both expressions equal, which is k = \ + 



i i i 



2 t 2(2/3+l) 2 ' 

leading to the claimed value. □ 

Lemma 13. There exists a path r' fromp to a point q' on the boundary of R via q, whose thickness 
is at most (1 + (2/3 + 1) 2 )|. 

Lemma 14. For any problem instance T>, there are at most crossing patterns 

between r' and p. 
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Proof. (Sketch) Let fi = (2/3 + l) 2 and v = \y e proceed as in the proof of Lemma 

Recall that previously we had 2 1,5 * x 2r • (3i x 3t x 4) 2r crossing patterns, but now we must use tne 
bounds that depend on (5 instead. What before was 2r now becomes fir, and the 1.5i terms now 
become vt (recall that 3t = 2 • 1.5i, hence they must be replaced as well). Making these changes in 
the previous expression, we obtain that the number of crossings is bounded by: 

2 vt x fir ■ (2vt x 2vt x 4)" r . 
Using that t < fir, and simplifying the expression, we get: 

< 2 u ^ r x fir ■ (A.vfir) 2 ^ r = 2^+log(Atr)+2^rlog(4!^r) 
Finally, we apply that both fi, v G 0(/3 2 ), and obtain the desired bound. □ 
Lemma 15. Any crossing pattern has at most 0(f3 r 2 ) forbidden pairs. 

Proof. This follows from the fact that r' and p traverse through at most 2fir cells of A(T>) (where 
fx = (2(3 + l) 2 ), hence at most such many crossings can occur. Each crossing generates a single 
boundary vertex, hence the number of forbidden pairs is bounded by 0((2fir) 2 ) = 0(/3 4 r 2 ). □ 

With these results in place, the rest of the algorithm remains unchanged: the only property of 
unit disks that is still used is the fact that they are connected, to be able to phrase the problem as 
a vertex cut in the region intersection graph. 

Theorem 4. Let T> be a collection of n connected (3-fat unit regions in M 2 , and let p and q be two 
points. Let r be a parameter. There exists an algorithm to test whether r(p, q) < r, and if so, to 
compute a path with that resilience, in 0(2-^ ,r )?i 3 ) time, where f(f3,r) G 0(/3 4 r 2 log(/3r)). 

Proof. As before, the running time is bounded by the product of the number of crossing patterns 



and the time needed to solve a single instance of the vertex multicut problem. By Lemmas 14 
and 
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these bounds now become 2°^ 2 ^°s(M) anc } 2 °^ r lo §(^)) , respectively. The product 
of both is dominated by the second term, hence the theorem is shown. □ 

5 (1 + e) -approximation 

An arrangement T> is said to have bounded ply A if no point p G M 2 is contained in more than 
A elements of T>. In this section, we present a polynomial-time approximation scheme (PTAS) for 
computing the resilience of an arrangement of disks of bounded ply A. The general idea of the 
algorithm is very simple: first, we compute all pairs of regions that can be reached by removing 
at most k disks, where we will set k = [4A/e 2 ]; the reason for choosing this value will become 
clear later. Then, we compute the shortest path in the dual graph of the arrangement of regions, 
augmented with extra edges for these pairs. We prove that the length of the resulting path is a 
(1 + ^-approximation of the resilience. As in the previous section, we first consider the case in 



which T> is a set of n unit disks in M 2 of ply A. In Section 5.4 we extend the result to /3-fat regions. 

Let A(V) be the arrangement induced by the regions of V, and let Gm%>\ be the dual graph of 
A(T>). Recall that Gm~d) has a vertex for every cell of A(T>), and a directed edge between all pairs 
of adjacent cells of cost 1 when entering a disk, and cost when leaving a disk. For any given k, let 
Gk be the graph obtained from G^(p) by adding, for each pair of cells A, B G A(V) with resilience 

at most k, a shortcut edge AL^ of cost r(A,B). 
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Figure 13: The point m is contained in at least c' + 1 disks of .4(2?). Therefore, one of the neighboring cells 
has ply at least c' + c + 1 

For a pair of cells of A(T>), we can test whether r(A, B) is smaller than k, and if it is, compute 
it, in 0(2^ k ^n 3 ) time (where f{k) = r 2 log r + o(r 2 log r)) by applying Theorem[3]to a point p 6 A 
and a point q £ B. Since there are 0(n 2 ) cells in *4(P), we can compute Gk by doing this 0(n 4 ) 
times, leading to a total running time of 0(2^ k ^n 7 ). Observe that this running time is polynomial 
in n, and exponential in k. In particular, it is a PTAS since k = 4A/e 2 . Again, we emphasize that 
the bounds are loose, and that our objective is to show the existence of a PTAS to the resilience 
problem. 

5.1 Analysis 

We will now show that the algorithm described above yields a (1 + e)-approximation. We begin 
with the following observation. 

Lemma 16. Let D 6 T>, where A(T>) has ply A, and let s,t be any two points inside D. Then the 
resilience between s and t inD is at most A. 

Proof. Let c be the number of disks containing either s or t (c > 1, since D contains both points). 
These c disks clearly must be removed. Now we analyze what other disks, not containing neither 
s nor t, may need to be removed too. For each other disk D\ (not containing both s and s) that 
needs to be removed in an optimal solution, there must be another disk D2 that intersects D\ and, 
together, separate s and t inside D. We call such a pair of disks a separating pair. 

Thus if the resilience is (c + c'), there must be at least c' disjoint separating pairs intersecting 
D. Moreover, since disks have unit-size, if two disks form a separating pair, at least one of them 
must intersect the center of D. Similarly, any disk containing s and s must contain the center of 
D as well. Fig. [13] illustrates the argument. 

Since the ply of A(D) is A, this implies that there can be at most A — c separating pairs, and 
thus the resilience is at most A. □ 

The previous lemma implies that in an optimal resilience path, if a disk appears twice, its two 
occurrences cannot be more than 2A apart (when counting the cells traversed by the path between 
the two occurrences of the disk) . 

To prove our result it will be convenient to focus at the sequence of disks encountered by a path 
when going from p to q. It turns out that such problem is essentially a string problem, where each 
symbol represents a disk encountered by the path. In that context, the thickness will be equivalent 
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Figure 14: (a) A string of 52 symbols, each appearing twice, (b) First, we identify a maximal 
set of segments bounded by equal symbols, and longer than A = 4. (c) Then, we segment the 
remaining pieces into segments of length k = 10. Red symbols are double-counted. 

to the number of symbols of the string (recall that we assume that p is not contained in any disk), 
and the resilience to the number of distinct symbols. 



5.2 Observations about strings 

Let S = (si . . . s n ) be a string of n symbols from some alphabet 21, such that no symbol appears 
more than twice. Let T be a substring of S. We define £(T) to be the length of T, and d(T) to 
be the number of distinct symbols in T. Clearly, \£{T) < d(T) < £(T). Let a and k be two fixed 
integers such that a < k. We define the cost of a substring T of S to be: 

a if T = (aXa) for some a G 21, string A such that a A, and £(T) > a 
ijj{T) = \ d(T) if £(T) < k 
^£{T) otherwise 

Note that, in the string context, d acts as the resilience, £ as the thickness, and tp is the 
approximation we compute. Intuitively, if T is short (i.e., length at most k) we can compute the 
exact value d(T). If T has a symbol whose two appearances are far away we will use a "shortcut" and 
pay a (i.e., for unit disk regions, by Lemma 16 we will take a = A). Otherwise, we will approximate 
d by I. Given a long string, we wish to subdivide S into a segmentation T, composed of m disjoint 
segments (i.e. substrings of S) T\, . . . , T m , that minimize the total cost ift(T) = Y^,i VK^i)- Clearly, 
4>(T) < £(S). 

Lemma 17. Let S be a sequence. There exists a segmentation T such that ip(T) < (1 + e)d(S), 
where e = 2^/ 'a jk. 

Proof. Let A be an integer such that a < A < k, of exact value to be specified later. First, we 
consider all pairs of equal symbols in S that are more than A apart. We would like to take all of 
these pairs as separate segments; however, we cannot take segments that are not disjoint. So, we 
greedily take the leftmost symbol s whose partner is more than A further to the right, and mark 
this as a segment. We recurse on the substring remaining to the right of the rightmost sffl Finally, 



we segment the remaining pieces greedily into pieces of length k. Fig. 5.2 illustrates the resulting 
segmentation. 

Now, we prove that the resulting segmentation has a cost of at most (1 + e)d(S). First, consider 
a symbol to be counted if it appears in only one short (blue) segment, and to be double- counted if it 
appears in two different short segments. Suppose s is double-counted. Then the distance between 
its two occurrences must be smaller than A, otherwise it would have formed a long (red) segment. 



m ln fact, we could choose any disjoint collection such that after their removal there are no more segments of this 
type longer than A. 
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(a) (b) 

Figure 15: (a) The optimal path ir* , achieving a resilience of 2. There is a segmentation of tt* of cost 
3, using the dashed shortcut, (b) A minimum cost path ir found by the algorithm. In this example, the 
resilience of tt is 3. 

Therefore, it must appear in two adjacent short segments. The leftmost of these two segments has 
length exactly k, but only A of these can have a partner in the next segment. So, at most a fraction 
X/k symbols are double-counted. 

Second, we need to analyze the cost of the long (red) segments. In the worst case, all symbols 
in the segment also appear in another place, where they were already counted. In this case, the 
true cost would be 0, and we pay a too much. However, we can assign this cost to the at least A 
symbols in the segment; since each symbol appears only twice they can be charged at most once. 
So, we charge at most a/A to each symbol. The total cost is then bounded by (1 + X/k + a/X)d(S). 
To optimize the approximation factor, we choose A such that X/k = cr/X; more precisely we take 
A = [V^cr] . Recall that we initially set k = \4a/e 2 ~\. □ 

5.3 Application to resilience approximation 

We now show that the shortest path between any p, q in Gk is a (1 + ^-approximation of their 
resilience. Let tt be a path from p to q in M 2 , and let S(tt) be the sequence that records every disk 
of T> we enter along tt, plus the disks that contain the start point of tt, added at the beginning of 
the sequence, in any order. Then we have | ^(vr) | = t(n). 

Lemma 18. For every path tt from p to q and every segmentation T of S(ir), there exists a path 
from p to q in Gk of cost at most ip(T). 

Proof. We describe how to construct a path in Gk based on T ■ For every segment T of T, we 
create a piece of path whose length in Gk is at most the cost of the segment ip(T). 

There are three types of segments. The first type are segments that start and end with the 
same symbol o, which corresponds to a disk D G T>. For those, we make a shortcut path that 



stays inside D, as per Lemma 16 The second type are segments whose length is at most k. For 
those, by definition, Gk contains a shortcut edge whose cost is exactly the resilience between the 
corresponding cells or A(T>). The third type are the remaining segments. For those, we simply use 
the piece of it that corresponds to T. □ 

Lemma 19. For any points p,q £ M 2 , the resilience of 3pQ k (p,q) is at most (1 + e)r{p, q). 
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Proof. Let tt* be a path from p to q of optimal resilience r* = r(n*) = r(p,q). Then, consider 



the sequence S(tt*), that is, the sequence of disks tt* enters. Now, by Lemma 17, there exists a 



segmentation T of S(tt*) of cost at most (1 + e)d(S(n*)) = (1 + e)r* . By Lemma 18, there exists 



a path in Gk of equal or smaller cost. Fig. |15(a) illustrates this. 



Now, consider the path tt our algorithm produces. The resilience of tt is smaller than the cost of 
tt in Gk, which is smaller than the cost of tt* in G^, which is smaller than 1 + e times the resilience 
of tt* . That is: 

r(vr) < cost Gfc (tt) < cost Gfe (tt*) < (1 + e)r(-TT*) = (1 + e)r* 

□ 

We conclude with our main result in this section. 

Theorem 5. Let T> be a set of unit disks of ply A in M 2 . We can compute a path tt between any 
two given points p,q £R 2 whose resilience is at most < (1 + e)r(p,q) in 0(2^^' £ ^n 7 ) time, where 
f(A,e) = 16^im + o(^im). 

Proof. The running time of the algorithm is dominated by the preprocessing stage: determining 
if the resilience between every pair of vertices of G^(p\is at most [~4A/e 2 ] . Since G^x)) has 0(n 2 ) 
cells, in total we execute the algorithm of Theorem at most 0(n 4 ) times, and we obtain the 
desired bound. □ 

5.4 Extension to Fat Regions 



As in Section 4.2 we now generalize the result to arbitrary /3-fat unit regions. For the purpose, we 



must extend Lemma H6l 

Lemma 20. Let D £ T>, where A(T>) has ply A, and let p,q be any two points inside D. Then the 
resilience between p and qinVis at most (2/3 + 1) 2 A. 

Proof. The resilience between p and q is upper-bounded by the number of regions that intersect D. 
We can give an upper bound using a simple packing argument. Since p and q belong to a /3- (unit) fat 
region D, they are both inside a circle C with center c and radius (5. Any other /3-fat region D' 
that interferes with the path from p to q must intersect C. Such an intersecting region, being also 
/3-fat, must contain a unit-disk whose center cannot be more than 2/3 away from c. Therefore all 
regions intersecting C have their unit-disks centered at distance at most j3 from c. Moreover, such 
disks are totally contained in a disk of radius 2/3 + 1 centered at c. As in the proof of Lemma 11 , we 
can show that at most (2/3 + l) 2 disjoint unit-disks fit into a disk of radius (2/3 + 1). Since the ply 
is at most A, the maximum number of unit-disks inside a disk of radius /3 in T> is (2/3 + 1) 2 A. □ 

As before, the rest of the arguments do not rely on the geometry of the regions anymore, and 
we can proceed as in the disk case. The only difference, is that the value a of doing a shortcut has 
increased to (2/3 + l) 2 A. 

Theorem 6. Let V be a set of unit disks of ply A m I 2 . We can compute a path tt between 
any two points p,q € M 2 whose resilience is at most < (1 + e)r(p, q) in 0(2-^ A '^ ,e )n 7 ) time, where 
feO(^flog(/3A/e)). 
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